This Page Is Inserted by IFW Operations 
and is not a part of the Official Record 

BEST AVAILABLE IMAGES 



Defective images within this document are accurate representations of 
the original documents submitted by the applicant. 

Defects in the images may include (but are not limited to): 



BLACK BORDERS 

TEXT CUT OFF AT TOP, BOTTOM OR SIDES 
FADED TEXT 
ILLEGIBLE TEXT 
SKEWED/SLANTED IMAGES 
COLORED PHOTOS 

BLACK OR VERY BLACK AND WHITE DARK PHOTOS 
GRAY SCALE DOCUMENTS 



IMAGES ARE BEST AVAILABLE COPY. 



As rescanning documents will not correct images, 
please do not report the images to the 
Image Problem Mailbox. 



per 



WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 




INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International Patent Classification 4 : 
H04N 1/46 



Al 



(11) International Publication Number: 
(43) International Publication Date 



WO 86/ 06907 

20 November 1986 (20.11.86) 



(21) International Application Number: PCT/US86/00823 

(22) International Filing Date: 23 April 1986 (23.04.86) 



(31) Priority Application Number 

(32) Priority Date: 

(33) Priority Country: 



730,627 
6 May 1985 (06.05.85) 
US 



(71) Applicant: EASTMAN KODAK COMPANY [US/US]; 

343 State Street, Rochester, NY 14650 (US). 

(72) Inventor: ALKOFER, James, Samuel ; 1668 Redman 

Road, Hamlin, NY 14464 (US). ' 

(74) Agent: CLOSE, Thomas, H.; 343 State Street, Roches- 
ter., NY 14650 (US). 



(81) Designated States: DE (European patent), FR (Euro- 
pean patent), GB (European patent), JP. 



Published 

With international search report. 



(54) Title: DIGITAL COLOR IMAGE PROCESSING EMPLOYING HISTOGRAM NORMALIZATION 



J4 r t2 



DATA a CONTROL BUS 




56\ 



DAC - J 




60^ 



SYSTEM 
CONTROL 



f38 



\kocx 

AVG. 
FILTER 



HS TOG RAM ^„ 
_, COMPILER ii 



JO 



40 



' CCN^RAST 
MEASUREMENT- 



FILTER 



it 
0" 
R £ 



46: 



Li 



H 50 



STATISTICS 
COMPUTER 



y2Q 



7? 



•44 



H 



36, 



INPUT 
CALIBRATION 
TABLE 



DIGITAL IMAGE 
PROCESSOR 



CONTRAST 
INTERVAL 
SELECTOR 



HISTOGRAM 
NORMALIZE 



CCNTRASJ 
ADJUST 



54 



152 



22) 



COLOR REPRODUCTION FUNCTION GENERATOR 



COLOR 
REPRODUCTION 
FUNCTION 
TABLE R 



OUTPUT 



(57) Abstract 



A digital color image processing method is characterized by generating color reproduction functions by normalizing 
samples of color values from the image, and processing the color components of the color digital image by applying the 
color reproduction functions to the respective color components of the digital color image. The method adjusts both tone- 
scale and color balance automatically. 
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DIGITAL COLOR IMAGE PROCESSING 
EMPLOYING HISTOGRAM NORMALIZATION 

Technical Field 

5 The invention relates to digital color 

image processing and more particularly, to digital 
color image processing wherein tone-scale' 
corrections are effected by histogram modification. 
Background Art 

10 In the practice of digital color image 

processing, an original color image, such as a 
photographic negative, is sampled periodically in 
three colors (e.g. green, red and blue) to produce a 
digital representation of the original color image, 

15 The digital color image is processed by applying 
digital image processing functions to improve such' 
image qualities as sharpness, tone-scale, and color 
balance. The processed digital color image is then" 
displayed on a display medium such as photographic 
• 20 film or paper. 

Fig. 2 is a schematic diagram illustrating . 
apparatus for digital image processing. Such 
apparatus includes an input device 1 for sampling 
the original color image in each of three colors, 

25 and analog- to-digital converters 2 for producing the 
digital color image in the three colors. Commonly 
employed input devices include drum and flat bed 
scanners, linear and area solid state image sensing 
arrays, and CRT and laser flying spot scanners, each 

^0 being provided with appropriate color filters to 
produce the color separations. 

Each digital color separation image is 
stored in a mass storage memory 3, such as a solid 
state memory frame-buffer, magnetic tape or disc 

35 storage device. A digital computet 4 applies the 
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various image processing functions to the digital 
color image to produce the processed digital color 
Image. 

.The digital computer 4 may comprise a main 

5 frame general purpose digital computer, or for 

higher spee)d operation, a digital computer specially: 
, configured for high speed digital processing of 
color images . 

The processed digital color image is 

10 converted to sampled analog form by 

digital-to-analog converters 5 and is displayed on 
an output device 6 such as a drum or flat bed 
graphic arts scanner, or a CRT or laser flying spot 
scanner. The elements of the. image reproduction 

15 apparatus communicate via a data and control bus 7. 

As noted above, among the processing 
functions performed by the digital computer are the 
improvement of the tone- scale and color balance of 
the color image. In the article entitled "Tone 

20 Correction of Color Picture by Histogram 

Modification" by Yoichi Miyake, Nippon Shashin 
Sakkaishi, V. 48(2), pp. 94-101, 1980, the author 
proposes a digital color image processing method 
wherein the tone-scale corrections are effected by 

25 modifying the histogram of color values of the green 
separation image. Color corrections are implemented 
by solving a conventional set of color masking 
equations of the form: . 

R f = auR + a l2 G = ai 3 B (1) 

30 G 1 = a 2 iR + a 2Z G = a 23 B (2) 

B 1 = a 3i R + a 32 G = a 3 3 B (3) 
where the matrix of color correction coefficients 
a^ are. determined primarily by the 
characteristics of the input and output media. 

35 . : ' / 
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In the field of optical photographic 
printing, it is well known to adjust the 
coefficients of the masking equations to 

correct for "problems' 1 in the color balance of the 

5 original, such as a color negative film balanced for 
daylight exposure that was exposed under Tungsten 
. lighting conditions. . The correction is usually 
based upon the hypothesis that the large area 
integral of the colors in any scene should add up to 

10 a shade near gray. Full correction, i..e. forcing 
the integrated colors in the processed image all the 
way to gray will correct for problems in the 
original image caused by such factors as film 
keeping, incorrect exposure, or an improper match 

15 between the illuminant and the film balance. On the 
other hand, for scenes that do not obey the 
assumption that all the colors, add up to gray 
(commonly called color subject failure scenes), full 
correction introduces unwanted color shifts in the 

20 printed image. Often some compromise between zero 
and full correction is employed in optical printers 
based upon the average characteristics of an 
expected population of photographic negatives. 

This approach results in a large percentage 

25 of acceptable prints, but has a shortcoming in that 
it may not produce an "optimum correction 11 for any 
: given original. Furthermore ,' there exists a class 
of "problem 11 originals that simply cannot be 
corrected by adjusting the coefficients of the 

30 masking equation. This class includes those 

originals in which the sensitometric curves Cross 
each other. This phenomena is caused, for example, 
by severe abuse to the original such as exposure of 
undeveloped photographic film to formaldehyde. 
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The objects of the present invention are 
directed to providing a digital color image 
processing method for adjusting the tone-scale and 
color balance of the processed image that overcome 
5 the shortcomings noted above. 
Disclosure of the Invention 

The objects are achieved according to the 
present invention by adjusting both tone-scale and 
color balance of a digital color image by modifying 
10 the histograms of each of three primary colors in 
the image . 

The method is based upon two principle 
observations regarding the statistical properties of 
the color values in a high quality color image. The 
15 first of these principles is that a truly random 

sample of color values (e.g. photographic density or 
log radiance) in a high quality color image will 
form a normal (Gaussian) distribution. The second 
principle is that the standard deviation of a random 
sample of color values is invar ient with 
wavelength. That is to say, a truly random sample 
of color values in one. color will have the same 
standard deviation as a truly random sample in 
another color. 

25 The first principle noted above implies 

that a function that normalizes a random sample of 
color values will serve well as a color reproduction 
function, assuming that any deviation from normality 
in the random sample was caused by some "problem 11 

30 with the original. The first principle combined 
with the second principle noted above implies - that 
color values in all three colors having an equal 
distance in their number of standard deviations from 
the. mean of their respective color distributions 

35 should always combine to produce a shade of gray. 
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Applyirig these principles to a digital 
color image processing method, a color reproduction 
function for each of three primary colors of the 
image is generated by normalizing a random sample of 
5 color values from each of the three colors. The 
color components of the digital color image are 
processed with the color reproduction functions for 
the respective colors to generate processed digital 
color image. 

10 xn a preferred embodiment of the invention, 

the sample of color values used to generate the 
color reproduction functions are taken from an image 
dependent "floating" contrast interval to optimize 
the randomness of the sample. The floating contrast 

15 interval is selected as a function of the shapes of 
the histograms of color values for the three primary 
colors in the contrast interval. This method of 
choosing a sample of image values is the subject of 
copending Patent Application Serial 

. 2U No. 730,627 • In the preferred embodiment 

of the invention, the color reproduction functions 
relate input color values to values of a standard 
normal variate. The term standard normal variate 
refers to a scale representing the number of 

25 standard deviations of a normal (Gaussian) 

distribution having a mean of zero and a standard 
deviation of one. Color components of the digital 
color image are processed by the color reproduction 
function to yield dimensionless Z values. The 

3° contrast of the processed image is determined by 
multiplying the Z values by a constant. The color 
balance of the. processed image is determined by 
adding respective constants representing the mean 
• density of the output medium for the respective 

35 colors, to the processed Z values; The basic method 
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of generating a tone reproduction function and 
determining contrast in this way is the subject of 
copending Patent Application Serial No . 730,630 * 

According to one embodiment, the function 

5 of histogram shapes employed to select the contrast 
interval is .the normality of the histogram in one of 
the three colors, weighted by the similarity of the 
shapes of the histograms in the other two colors to 
the one. Iri another embodiment, the function of 

1° histogram shape is the similarity in shapes of * all 
the color histograms in a given contrast interval. 
' z In one embodiment of the invention, the 

multiplicative constant employed to adjust the 
tone-scale of the processed image is generated as a 

15. function of the standard deviation of the sample of 
color values in one of the colors (i.e. green). 

In an embodiment where the digital image 
was generated by sampling photographic film having a 
film grain structure, the image is sampled at full 

^ti system resolution in three primary colors, then the 
color values in all three colors are block averaged 
prior to generating the color reproduction functions 
to reduce the effects of grain on the sampling 
process. After the color reproduction function are 

25 generated, the full resolution digital color image 
is processed. 

As a further refinement of the method, the 
total number of color values in each contrast 
interval is counted, and only the contrast intervals 

30 having a number of values greater than some 

predetermined minimum number are considered for 
generating the color reproduction functions. This 
is done to insure a statistically significant sample 
of color values. 

35 / 
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Brief Description of the Drawings 

. The invention is described with reference 

V 

to the drawings, wherein: 

Fig* 1 is a schematic diagram, illustrating 
5 a scanning color photographic printer for 

carrying out the digital color image processing 
method according to the present invention; 

Fig. 2 is a schematic diagram illustrating 
generic prior art digital color image processing 
10 apparatus useful for practicing the present 

invention; 

Fig. 3 is a schematic diagram showing the 
organization of the histogram memories; 

Fig. 4 is a graph illustrating a color 
15 value histogram from one of the contrast 

intervals shown in Fig. 3; 

Figs. 5-8 are flow charts illustrating the 
method of generating the color reproduction 
functions according to the present invention; and 
20 Fig. 9 is a graph showing the form of the 

color reproduction functions generated according 
to the method outlined in Fig. 8. 
Modes of Carrying Out the Inventions 

Turning now to Fig. 1, an example of a 
25 scanning color photographic printer used to practice 
the present invention will be described. 

The input device includes a light source 10 
and lens 12 for projecting an image of color 
negative film 14 onto three solid state image 
30 sensing arrays 16, such as CCD image sensors. The 
image is directed to the image sensing array 16 by 
dichroic beam splitters 18 through color trimming 
filters 20, to form a red, green and blue separation 
image on the respective image sensors. A film 
35 advance motor 21 advances the film in the printer. 
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The image sensors 16 scan the color photographic 
negatives to produce three color separation signals 
R, G and B. The signals thus produced are supplied 
to analog-to-digital converters 22 that each produce 

5 an 8-bit output code representing one of 256 

possible sigWl levels for each sample point in the 
three colors .• 

The digital color image signals are stored 
in a digital image storage memory 24, such as a 

10 magnetic tape, disc, or solid state semiconductor 
memory. The digital color image signals are 
processed in a digital computer 26, a DEC 2060 
mainframe computer was used. 

The digital computer 26 is programmed to 

15 include a color reproduction function generator 28 
and a digital image processor 30 that applies color 
reproduction functions generated by color 
reproduction function generator 28 to the digital 
color image. The color reproduction function 

20 generator 28 receives the digital color image from 
the digital image storage device 24 via data and 
control lines 31 and generates a color reproduction 
function table 32 for each color. The color 
reproduction function generator 28 also generates a 

25 multiplicative constant a z for adjusting the 

contrast of the processed image as described below 
and supplies the constant to three output 
calibration tables 34. 

Prior to storing the digital color image in 

30 digital image store 24, each input signal level is 
converted to an input color value, such as 
photographic density, by a known scanner calibration 
function. The calibration function for each color 
is implemented in the form of a look up table 36 in 

35 the digital image processor 30, to convert each 
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8-bit input value to an 8-bit value representing the 
color negative density at the scanning point. 

The color reproduction function generator 
28 includes digital filters 38 for performing a 
5 block average of the color values of the color 
components of the digital color image. Digital 
filters 40 measure the contrast of the image around 
each block averaged color value in each of the three 
colors. A histogram compiler 42 compiles the block 

10 averaged color values from a plurality of contrast 
intervals in each color, in a plurality of histogram 
memories 44. 

Fig. 3 shows, in a graphic way, the 
organization of each of the histogram memories 44. 

15 There are twenty contrast intervals for each color, 
having a width of .04 log contrast units each. The • 
width of the contrast intervals is chosen to be 
approximately twice the. minimum visual* log contrast 
threshold. The width of the contrast intervals 

20 represents a tradeoff between randomness of sampling 
(the narrower the interval, the greater the 
randomness) and achieving a sample size that is 
statistically significant. The 256 color values are 
divided into 80 color value (density) bins, for a. 

25 resolution of .05 density units per density bin. 
Counts are accumulated in the appropriate density 
bins in the histogram memories until all of the 
color values from the digital color image are 
counted. Fig. 4 shows a graphic example of one of 

30 the color histograms representing one of the 
contrast intervals in one of the three colors. 

Returning to Fig. \ y a statistics computer 
46 in the color reproduction function generator 28 
computes the statistical moments y of the 

35 -histograms H in the histogram memories 44 as 
described below. 
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A contrast interval selector 50 receives 
the statistical moments p from the statistics 
computer 46 and selects a contrast interval on the 
basis of the statistical moments of the color value 

5 distributions in the contrast intervals and supplies 
the histograms H from the selected contrast interval 
to a histogram normalizer 52. The histogram 
normalizer 52 normalizes the histograms from the 
selected contrast interval to generate the color 

10 reproduction function lookup tables 32 that are used 
by the digital image processor 30. 

A contrast-adjustment cdmputer 54 receives 
the standard deviation a Q of the green color 
values in the selected contrast interval and 

15 generates a multiplicative constant used to 
determine the contrast of the processed image in 
. output device calibration tables 34. All the 
processed color values in all three colors are 
multiplied by this same constant. 

20 The color reproduction function lookup 

tables 32 relate each of the 256 possible input 
values in the respective colors to one of 256 
possible output values. After the color 
reproduction function lookup tables 32 have been 
* 25 generated, the digital image processor 30 processes 
the digital image by applying the color reproduction 
functions to the respective color components of the 
digital color image. The output calibration 
functions are then applied to the respective color 

20 components of the processed digital color image. 

The processed digital color imagers converted back 
to analog form by digital-to-analog converters 56. 
The processed analog signal is then applied to the 
output scanning device 58 to reproduce the image. 

35 ' ~ / 
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The output device 58 is a drum-type scanner having 
red, green and blue laser light sources that are 
modulated by the respective analog color signals to 
expose a light sensitive medium such as color 
photographic paper. 

Syltem control electronics 60 controls and 
coordinates the operation of the input and output 
scanning devices and the digital image processing 
computer 26 . 

The method of generating the color 
reproduction functions will now be described in more 
detail with reference to the flow charts of Fig. 5-7, 

Referring first to the flow chart of Fig. 
5, the processing steps performed on the color 
15 digital image to generate the respective color 

reproduction functions will be described. All three 
color components are processed through these steps. 
First a block average of the color values of the 
sampled digital image is formed. This is 
accomplished by applying a digital filter to the 
digital image color values of the form: 



,20 



25 
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(4) 



This .averaging is performed by the block average 
filters 38 shown in Fig. 1, to reduce the effects of 
30 film grain and other sources of noise on the color 
value statistics. 

Next, a digital filter representing a 
Laplacian operator of the form: 



35 
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(5) ' 



• .? 

is applied to the' block averaged color values of the 
digital image to measure the contrast of the image 
at each block averaged sample point. This operation 

j^q is performed by the contrast measurement filters 40 
..in Fig, 1. The Laplacian operator. has the 
properties of exhibiting no response. to uniform 
areas or linear gradients, and exhibits a response 
only to changes in gradients, the Laplacian^ 

15 operator works well in measuring the contrast of the 
- modulated color components of the image, however it 
is to be understood that other contrast measuring 
filters may be employed with the present invention. 

The respective histograms in each contrast 

20 interval for each color component are. compiled as 
discussed above , and their statistics computed. A 
contrast interval is selected on the basis of the 
shape of the histograms of the respective color 
values in the interval and the histograms of color 

25 values in the selected contrast interval are 
normalized to generate the respective color 
reproduction functions 

Fig. 6 is a flow chart showing the steps 
involved in compiling the histogram statistics. The 

30 raw. moments taken about the mean, are 
computed as follows: 

N 

y k = I ( x L - x ) k } (3) 

N. i=l 

35 
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where 

N is total number of samples in the 
contrast interval; 

x^ is a color value; and 
5 x is the mean color value. 

The standardized central moments y' k are 
calculated as follows: 

= , where a = / pa (7) 

10 a k 

The coefficient of symmetry (coefficient of 

skewness) for each distribution is then represented 
as 

81 - (u's) 2 (8) 

-L5 and the coefficient of peakedness (kurtosis plus 3.) is 
represented as 

82 = \i\ (9) 
. Referring now to Fig. 7A, the steps 

involved in selecting the contrast interval 
according to one embodiment of the invention will be 
described. One selection criteria that has been 
found to work well is a combination of the normality 
of the sample of green color values and the sameness 
of the distribution of the red and blue color 
25 values. A function called RANK is defined as 
follows: 

(10) 

RANK = /B 2 i G + (8 2 Q -3) 2 + | B i R - B i Q I + I B 2 R - B zq I 



20 



30 



35 



+ I Bi B - Bi Q I + I Bz B - Bag I 

where 8 in is the skewness of the distribution of 
green color values in the contrast interval being 
evaluated, etc. 

The contrast interval with the lowest value of RANK 
is selected. 
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In a presently preferred, alternative 
embodiment of the invention, (shown in Fig. 7B) the ■ 
contrast interval selection criteria involves the 
similarity of shape of the three color value 
5 distributions in the contrast interval. The 

selection criteria compares the skewness of the 
color value distributions in each of the contrast 
intervals as follows : 

(ID 

10 SUMi(N) = | Bi R - Bi G I + I Bi B - 6i 6 I + I Bi R .- B i B I 

. where SUM i (N) is a measure of the difference in 
skewness between the three color value distributions 
in the Nth contrast interval. The contrast interval 
having the most similar distribution, i.e. the 
15 lowest value of SUMi (N) is identified. 

Next, the peakedness of the three color 
value distributions in each contrast interval is 
compared as follows: . 

(12) 

20 SUM 2 (N) = | B* R - f*2 G | + I Bz B * Ba Q I + I B i R - Big I 

where SUM 2 (N) is a measure of difference in 
peakedness between the three color value 
distribution in the Nth contrast interval. The 
contrast interval having the most similar 

25 distributions, i.e. the lowest value of SUM 2 (N) 
is identified. 

Finally, the lowest contrast interval (i.e. 
the interval representing the lowest contrast) 
between the two identified contrast intervals is 

30 selected. 

Optionally, to insure that there is a 
statistically significant number of color values in 
the contrast interval, a check on the total count of 
color values in the interval is performed. If the 
35 total count in each color is less 7 than some 
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predetermined number, say 1000, the contrast 
interval is not selected. This check is continued 
until a contrast interval having at least the 
required minimum number of color values is chosen. 

5 Turning now to Fig. 8, the steps involved 

in normalizing the color histograms in the selected 
contrast interval to generate the color reproduction 
funcations will be described. To normalize a 
histogram, the standard normal variate Z for all 80 

10 color value bins in the contrast interval is 
computed. First however, an average smoothing 
operation is performed on the selected histogram to 
remove any spikes. The smoothing is performed on 
the counts in 3 consecutive color value bins as 

15 follows: 

h. « 1/3 (h'. m l + h^ + h'. + x ) (13) 
where h'^ is the count in bin i and 

h. is the smoothed value. 

Next, the standard normal variate Z is 
^0 calculated for the smoothed values of the histogram 
as follows (from Approximations for Digital 
Computers , Hastings C, Princeton Univ. Press.): 



ao + a i t 



z. - t. 



j (14) 



25 " 1 + bitj + bztj 2 

where 



t . = / In (1/Pj 2 ) 

a 0 = 2.30753 b t = .99229 

30 ai = 0.27061 b 2 =. .04481 

The probability Pj for each of Che 80 bins is 
given by 

35 P. = i=1 " , (15) 



80 

Z h. 
i=l 
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where 

h L are the smoothed counts in the ith 
color value bin, and 
j = 1 - 80. 

5 Next, the Z values are linearly 

interpolated rrom 80 to 256 values to provide a Z 
value for each of the 256 possible scanner inputs 
represented by the 8-bit digital code. Finally the 
.256 Z values are stored in the color reproduction 
-10 function lookup tables 32. This process is applied 
to each of the three samples of color values in the 
selected contrast interval. 

Fig. 9 is a graph showing the form of the 
color reproduction functions produced by normalizing 

15 the samples of color values. In the lower right 

quadrant of the graph, a curve labeled 70 represents 
a standard normal distribution showing the 
probability of the occurrence of a value plotted 
against the standard normal variate Z. In the upper 

20 left quadrant of the graph, the curve labeled 72 
represents the sample of color values from the 
image, plotted against relative probability of 
occurrence. The central ordinate of the graph 
relates the relative probability Pj of the color 

25 value distribution to Z values according to the 
relationship defined by equation 14. The color 
reproduction curve, labeled 74, maps the Z values on 
the ordinate to the same Z values on the abcissa. A 
color value scale on the far upper right of the 

30 diagram, congruent to the color value scale on the 
far left, shows how the color reproduction function 
relates color values to Z values. 

After the three color reproduction function 
lookup tables are generated, all of the color values 

35 of the image are processed by applying the 
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respective color reproduction functions to them. At 
this point, the processed color samples from the 
image are dimensionless quantities representing the 
Z values. 

5 To recover the processed image, these 

dimensionless quantities are given magnitude with 
respect to the original image and the output medium 
by multiplying the Z values with a multiplier 
that adjusts the contrast of the processed image. To 

10 adjust the color balance of the image, a constant 
term is added to each of the three color values. 
The constant term relates the Z values in each color 
to the mean density of the output medium for the 
respective color, thereby causing equal Z values in 

!5 all three colors to be reproduced as a shade of gray. 
Appropriate values, for the multiplier 

a that adjusts contrast and the additive 

. z . • 

■ constants that determined- the color balance of the 
processed image are determined as follows. The 

20 multiplier is computed based on the statistics of 
the green color values alone, /but is applied to all 
three colors. The intrinsic contrast of natural 
scenes can be quantified in terms of the standard 
deviation of log reflectance of edges in the scene 

25 or the density representation of those log 

reflectances in the photographic image. On the 
• average the approximate relationship between the two 
is given by: 



where: 



G = average gradient of the photographic 

film (relates a R to some specific 
reproduction medium contrast)' 



35 
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standard deviation of log reflectance 
based on a large number of original 
scenes 

standard deviation of density 

!> Typical values for c R and G for color 

photographic film are .31 and .68 respectively, such 

that o D is .21. Departures from this average 

contrast must be compensated. A general equation 

may be stated as:-' 

10 - a = m . f(a q ) +b (17) 

z » 

where: 

a = individual scene standard deviation, 

from the selected contrast interval' 

m and b are system dependent constants and 

15 £( a ) is -some function of the sample standard 

s 

deviation 

a = the multiplier applied to the values 

obtained from the color reproduction 
function. 

A simple and satisfactory implementation is 
obtained from: 

b - a D . (1.0 - m) (18) 
a z = a . a s 4- b (19) 

where : 

m is typically between .6 and .8. 

The sign of a z is positive if the 
reproduction. has the same polarity as the original 
image (negative- to-negative or 

positive-to-positive). If the reproduction has a 
polarity of an opposite sense with respect to the - 
original, e.g., negative-to-positive, then the sign 
of a z is negative. 

Note that the adjustment of the contrast in 
(19) does not affect scenes having ; the average 
contrast (a D ) , nor does .it affect the mean value 
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o£ tone in the scene, since in terms of Z (the 
normal variate) the average value remains 0. This, 
is not only a computational convenience in terms of 
adjusting the contrast, but also in adjusting the 

5 color balance of the processed image with respect to 
the reproduction medium. For example, if a color 
negative image is to be printed directly onto color 
photographic paper, the log exposure for the desired 
mean paper density for each color is simply added to 

10 the translated, contrast adjusted values of the 
respective colors. The complete calculation is 
given by: " . 

lo § E ZD " ~ a z ' Z D + l0g E A (20) 
where : 

15 log E A - log exposure required to obtain 

the aim paper density 
= translated Z value for some input 

density in the original image 
log E ZD = log exposure for Z Q . 
Industrial Applicability and Advantages 

The digital color image reproduction method 
according to the present invention is useful in the 
graphic arts and photographic printing fields to 
provide automatic tone-scale adjustment and color 
25 balance of digitally processed images. The method 
is advantageous in that a greater percentage of high 
quality images are produced automatically, without 
the need for operator intervention, than by the 
methods of the prior art. 
30 By processing the color values according to 

the present invention to yield dimensionless 
quantities relating to the statistical distribution 
of color values, the tone-scale and. color balance of 
the processed image are easily adjusted by the 
35 simple acts of multiplying the processed values by a 
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constant to adjust the tone- scale, and adding a 
constant to the tone values in each color 
representing the mean density of the output medium 
to adjust the color balance. Because the processing 
5 o£ the original digital color image does not involve 
dimensioned quantities until the last step, the 
processing is completely independent of the output 
medium until the last step is performed. The method, 
therefore has the advantage that changing from one 
10- output medium to another is easily accomplished by 
I simply changing the constants representing that mean 
t density of the output medium desired. Under and 
r. over exposure are automatically compensated, since 
• the correction depends on the expectation of 
15 normally distributed density, and not on a priori 
information about film sensitometric 
characteristics. An additional advantage of the 
method according to the present invention is that 
color corrections are automatically performed by the 
20 method that would be impossible using the prior art 
methods . 
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Claims 

1. A digital color image processing 
method characterized by: generating color 
reproduction functions by normalizing samples of 
5 color values from the image, and processing the 
color components of the color digital image by 
applying the color reproduction functions to the 
respective color components of the digital color 
image . 

10 2. The digital color image processing 

method claimed in Claim 1, wherein said samples of 
color values from the image are selected from an 
image dependent contrast interval based upon the 
similarity in shape of the distribution of color 

15. values in the contrast interval. 

3. The digital color image processing 
method claimed in Claim 1, wherein said color 
reproduction functions relate input. color values to 
values of a standard normal variate Z, and wherein 

2^ the tone-scale of the processed image is determined 
by multiplying the Z values by a constant, and the 
calibration is accomplished with respect to a 
particular output medium by adding a constant to 
processed color values representing the mean density 

25 of the output medium for the respective color. 

4. The digital color image processing 
method claimed in Claim 3, wherein said contrast 
determining constant is a function of the standard 
deviation of one of the samples of color values. 

30 5. xhe digital color image processing 

method claimed in Claim 4, wherein said function is 
of the form: 

a = m • a + b 
z s 

35 



BNSDOCID: <WO 8606907A 1 _L> 



WO 86/06907 



PCT/US86/00823 



-22- 

when 

a z is the contrast adjusting constant, 

a is the standard deviation of the 
g 

sample of green color values, and 
5 m and b are system dependent constants. 

6 J A digital color image processing 
method for adjusting the tone-scale and color 
balance of the processed image, characterized by: 

a) generating a color reproduction 

10 function for each of three primary colors by 

normalizing a. random sample of color values from 
' the image in each of the three ' colors ; 

b) processing the color values of the 
image in the three primary colors using the 

15 respective color reproduction functions for the 

colors to generate processed values representing 
* a number of standard deviations from the mean 
for each color value; 

c) multiplying the processed values by a 
^0 constant to adjust the tone-scale of the 

processed image; and 

d) adding a constant representing the 
mean density of the output medium for each of : 
the three colors to the respective color values 

25 of the three colors ' to calibrate the output 

image with respect to a particular output medium, 

7 . The invention claimed in Claim 6 , 
wherein the samples of color- values employed to 
generate the color reproduction functions are from 

30 an image dependent "floating 11 contrast interval. 

8. The invention claimed in Claim 7, 
wherein the "floating" contrast interval is selected 
according to predetermined criteria relating to the 
shapes of the histograms of the color values in the 

35 contrast interval. 
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9. The invention claimed in Claim 8, 
wherein the predetermined criteria is a function of 
the normality of the histogram of one of the colors, 
and the similarity of the shape of the histograms of 
the other two colors to said one. 

10j The invention claimed in Claim 8, 
wherein the predetermined criteria is a function of 
the similarity of shape of the histograms of the 
three color values in the contrast interval. 

11. The invention claimed in Claim 6, 
wherein the constant multiplier employed to 
determine the tone-scale of the processed image is 
generated as a function of the standard deviation of 
the sample of color values of one of the colors. 

12. The invention claimed in Claim 7, 
wherein the contrast interval has at least a minimum 
number of color values. 

13. Apparatus for processing a digital 
color image to adjust the tone-scale and color 
balance of the processed image comprising: 

a) means for scanning the image in three 
colors to produce color values and contrast 
values in each color; 

b) means for selecting a random sample of 
color values from each color; 

c) means for generating a color 
reproduction function for each of the three 
colors by normalizing the random sample of color 
values in each color; 

d) means for processing the color values 
of the image in the three primary colors using 
the respective color reproduction functions for 
the colors, to generate processed values 
representing a number of standard deviations 
from the mean for each color value; 
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e) means for generating a contrast 
adjusting constant based on the standard 
deviation of the sample of one of the colors; 

f) means for multiplying the processed 
5 values in all three colors by the contrast 

adjustir/g constant to adjust the tone-scale of 
the processed image; 

g) means for adding a constant 
representing the mean density of the output 

^0 medium for each of the three colors to the 

respective color values of the three colors to 

7 calibrate the processed image with respect to a 

particular output medium. 

14. The invention claimed in Claim 13, 
15 wherein the means for selecting a random sample of 

color values from the image selects said sample from 
an image-dependent "floating 11 contrast interval. 

15. The invention claimed in Claim 14, 
wherein said means for selecting a random sample 

20. from a "floating 11 contrast interval selects said 
contrast interval according to predetermined 
criteria relating to the shapes of the histograms of 
the color values in the contrast interval. 
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